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Abstract 

We examine the hydro dynamic phonon spectrum of a monolayer film adsorbed on the wall of 
a spherical pore. Due to the boundary conditions, the monolayer film exhibits a discrete phonon 
spectrum. The corresponding density of states per unit frequency is thus a set of delta functions 
and the heat capacity exhibits Arrhenius behavior at low temperatures. At high temperatures, the 
heat capacity approaches the C ~ T 2 behavior of a two-dimensional monolayer film. Our results for 
the spherical surface film are compared to previous calculations for films confined to a cylindrical 
surface. 
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Monolayer films have been extensively studied, both theoretically and experimentally, in 
part because they exhibit reduced dimensionality and novel phases [|I|, 0]. For example, low 
pressure adsorption on graphite usually yields a two-dimensional (2D) monolayer, whereas 
adsorption in carbon nanotubes can result in quasi-one-dimensional behavior || |], ^, 0. 
As a result of the reduced dimensionality, the phonon spectrum and the thermodynamic 
properties of such films are altered with respect to bulk properties. The basic reason is 
that the substrate imposes its boundary conditions on the adsorbate. For example, previous 
calculations of Vidales, Crespi and Cole (VCC henceforth) show that the heat capacity of a 
monolayer film adsorbed on the inner wall of a carbon nanotube exhibits a crossover from 
quasi-one dimensional to two-dimensional behavior when the thermal phonon wavelength 
A = hs/3 approaches 27iR, the circumference of the cylindrical adsorbed film (s is the sound 
speed and /3" 1 = ksT) 0. Due to the cylindrical geometry, finite in the radial direction 
and (assumed) infinite in the axial direction, the surface-parallel phonons have a dispersion 
relation: 

oo nk = syjk* + {n/Rf (1) 

where k is the quasi-continuous wave vector (parallel to the cylinder axis) and n is the 
azimuthal quantum number. 

Here we consider the analogous case of a thin film adsorbed on the inner surface of a 
spherical cavity. The outline of this paper is as follows. We derive the in-plane phonon 
dispersion relation for the spherical surface film and compare the density of states with 
corresponding results for the cylindrical and planar cases. This leads us to anticipate some 
of the heat capacity results. Finally, we identify a crossover behavior in the heat capacity 
and compare this with some experimental measurements. 

In the derivation of the phonon dispersion relation, we adopt a method similar to that 
of the VCC paper concerned with cylindrical films J7|. It is noted first that at low T, 
the surface-normal vibrations can be ignored since they acquire a gap at roughly the first 
surface-normal excitation energy of the isolated atom (~ 60 K for 4 He on graphite). Thus, 
in a liquid monolayer film, only the surface-parallel phonons are excited at low T. These 
are described by a continuous two-dimensional wave vector k for a planar film and a one- 
dimensional wave vector in the cylindrical case. Since a spherical film is confined (effectively 
zero dimensional), we expect the quantum number equivalent to a wave vector k to have 
discrete values. The derivation of the spectrum follows from the scalar Helmholtz equation, 

2 



describing the velocity potential \1>: 

V 2 ^ + k 2 ^ = (2) 

where k = uj/s. This yields solutions of the form \l/(cos#, 0) = P / m (cos6 l )$(0), where $(0) = 
e'""^, m — 0, ±1, ±2, ... and are the associated Legendre polynomials for integers I > \m\. 
The boundary condition imposed on the Legendre polynomials (to be finite when cos 9 = ±1) 



yields the expected discreteness of the eigenvalue k = yl(l + l)/R. The dispersion relation 
for the spherical film follows as: 



ui = sJUJ + l){s/R) (3) 

The corresponding density of states per unit frequency g(oj) is thus given by a sum of delta 
functions with multiplicity given by the m-degeneracy (21 + 1) of each energy level: 

9{») = E^V) = E( 2 ' + !)% - <*) ( 4 ) 
i i 

We want to compare the spherical film results with the cylindrical and planar densities 
of states . For the cylindrical film, the number of states within the interval did is: 

g cy \u)duj = 2^dk (5) 

where Lj (2tt) is the number of states in the ID reciprocal space range dk, and 2 arises from 
the (±&) degeneracy. Using the dispersion relation (Eqn.l), the cylindrical density of states 
corresponding to a specific band (n) is: 

£V) = (6) 
™ yfl - (sn/uR) 2 

where 0(x) is the Heaviside step function (1 for x > and for x < 0). The total density 
of states is the sum over all bands: 

oo 

g cyl (") = E afiuj) (7) 

n=— oo 

At low T, the azimuthal excitations (n > 1) are frozen out due to the gap E\o = hs/R (see 
reference @). Thus, the dispersion relation (Eqn.l) becomes u = sk when n=0 and the 
cylindrical density of states becomes constant go(ui) = — . The planar film has also a linear 
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dispersion relation u = sk, but k is then a 2D wave vector. The corresponding density of 
states is: 



9 2D (") = ^3 (8) 



A uj 

where A is the area of the surface. 

The spherical, cylindrical and planar densities of states are displayed in fig. 1. In order to 
get dimensionless numbers, we plot the density of states relative to the planar case for the 
same area (A = 2-kRL in the cylindrical case and A = AttR 2 in the spherical case). Notice 
that the largest density of states at low uj corresponds to the cylindrical film, followed by 
the planar film and finally, by the spherical film. This implies that the cylindrical film will 
exhibit the largest heat capacity (per unit area) at low T, followed by the planar and then, 
the spherical film. At high uj, however, the number of states in an interval duj is the same 
for all types of films considered, as shown in the appendix. Therefore, the heat capacities of 
the spherical, planar and cylindrical films should approach asymptotically the same values 
at high T. We note that the "high T" limit in this elastic regime differs from the very high 
T limit in an experiment. That is, the law of Dulong and Petit corresponds to a very high 
T constant behavior of C. The crossover to the latter regime occurs when kT exceeds the 
Debye temperature. We do not consider this regime of T here. 

We turn our attention to the thermodynamic properties of the spherical film. Using Eq.3, 
the energy of a phonon is: 

Ei = (9) 

There is no mode for / = 0, and thus the spherical film possesses a gap E x = (hsy/2)/R. 
The total energy of the phonon system is: 

F - V {21 + l)El (10) 
1=1,2... ex P\PEi) - 1 

where 21+1 is the degeneracy of the energy levels. In the high T limit, the sum can be 
replaced by an integral, which can be evaluated analytically (because 1 <C I for the predom- 
inant modes): 

hs roo 2fdl AC(3) 

sph ~ R Jo exp([3hsl/R) - 1 vr/^V 1 1 

where A = AttR 2 is the film area and ( is the Riemann zeta function. As in the case of the 
cylindrical film, the high T heat capacity per unit area of the spherical film is the same as 
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that of a planar film: 

C2d(T) 3C(3) 

k B A TT^flh 2 1 ' 

For simplicity, we define a reduced temperature: 

t = R/((3hs) (13) 
and the dimensionless heat capacity takes the form: 

As in the case of the cylindrical film JF) , the heat capacity is seen to be a universal function of 
the reduced temperature (the right side of Eq.13). Fig.2 shows the dependence on T of the 
spherical and cylindrical heat capacities, compared to the planar heat capacity (C '/ C p i anar ) . 
Notice that the planar film is the high T limit for both the spherical and cylindrical films, 
as expected. However, the heat capacity of the cylindrical film approaches it from above, 
whereas the heat capacity of the spherical film approaches it from below. This behavior 
was expected from the previous discussion of the densities of states. We note that the 
asymptotic approach to the planar limit occurs at a lower T in the cylindrical case than 
in the spherical case. This can be explained in at least two ways. One is that a cylinder 
(having one infinite radius of curvature) is closer to a plane than a sphere is. The other way 
is that the density of states (Fig.l) of the cylinder is closer to that of a planar surface than 
is the density of states for the sphere. In analogy to the crossover behavior of the cylindrical 
film 0, one may speak of a crossover temperature (T CT0SS ) from a zero- dimensional system 
to a two-dimensional system, defined as that temperature where the heat capacity is about 
half of that of a planar film. In the present case, from fig.2, kBT cross ~ At low T, the 
spherical heat capacity has an activated behavior: 

C sph /k B = 6r 2 exp(-V2/t) (15) 

We do not know of any experiments related to adsorption in spherical pores, but the irreg- 
ular confined pores of a disordered material may conceivably be approximated as spherical 
pores. Crowell and Reppy have investigated A He films adsorbed in two porous glasses, Vycor 
and aerogel, using high-precision torsional oscillator and calorimetry techniques 0, |9|, |10|| . 
At some critical coverage, the 4 He film undergoes a transition from nonsuperfluid to super- 
fluid behavior. In the nonsuperfluid regime, the heat capacity exhibits a linear dependence 



at large T, but C/T drops rapidly at low T, so a crossover temperature between these two 
behaviors is defined and calculated to be about 100 mK for the lowest coverage. In the con- 
text of our calculations, a working hypothesis is that the pore geometry is like an American 
football, i.e. an ellipse of revolution. This could show the behavior of a cylinder (C ~ T) 
at modest T, but exhibit a gap at low T, as in the case of a spherical geometry. We may 
estimate the heat capacity at low T with our spherical pore model, and use the cylindrical 
pore for higher T. For a pore diameter of about 70 A, as in the case of Vycor pores, we 
assume that the A He atoms form a spherical film of about 60 A diameter, so there are about 
10 3 He atoms in a pore. Then the heat capacity per mole is of the order of mJ/K, the 
same order of magnitude as the experimental heat capacity of the nonsuperfluid A He film. 
However, the calculated heat capacity exhibits the characteristic activated behavior seen in 
-Eg. 15 and does not follow the general trend of the experimental measurements, which looks 
more like a power law behavior. Also, the heat capacity experimental values at higher T 
can be obtained with our cylindrical pore model but only with an unreasonably large speed 
of sound. However, the crossover temperature given by our spherical model is of the order 
of 150 mK, somewhat comparable to the experimental crossover of 100 mK. Thus, it may 
happen that the geometry of the Vycor pores produces boundary conditions for the film, 
whose effect changes with increasing temperature, leading to a crossover behavior. 

In conclusion, the substrate geometry imposes constraints on the phonons of monolayer 
films adsorbed within a spherical pore. The first effect of the confinement is the discreteness 
of the phonon spectrum. Because of a low ui gap in the spectrum, the specific heat exhibits 
an essential singularity at T=0, whereas at high T, it asymptotes to that of a planar film 
(C ~ T 2 ). A somewhat different crossover behavior occurs in the case of a cylindrical film, 
which has quasi-one-dimensional behavior at low T (C ~ T) and planar (2D) behavior at 
high T. This pair of solutions may explain some of the heat capacity features of monolayers 
adsorbed in disordered porous materials, whose various pores have different geometries and 
sizes. 

We are grateful to the National Science Foundation, which has supported this research. 
We thank Paul Crowell and Moses Chan for helpful communications. 
APPENDIX 

This paper states that the density of states at high ui is the same for all three problems 
discussed here. We will give now a proof of this statement. Let us start with the case 
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of a spherical film, and choose an interval du between u\ = \Jl(l + l)(s/R) and u 2 = 



(I + I) (I + 2)(s/R), corresponding to consecutive values of 1 when I ^> 1. The number of 
states in this interval is: 



/ g sph (u)du = {21 + 1) / 5(u -u t ) = 2l + l (16) 



Thus, in the limit of high u, and correspondingly high 1, the number of states in an interval 
du is of the order 21. 

For the planar film, we consider the same interval du, and the number of states in this 
interval is: 

fs?^ = 2(,+1) (17) 

where the area of the planar film equals the area of the spherical film A = AnR 2 . For high 
1, the number of states of the planar film in this interval du is therefore of the order of 21, 
equal to that of the spherical film. 

The last case is that of a cylindrical film. The number of states corresponding to the 
interval du is: 

r^ 2 L ' <d(u-sn/R) L ' r<^ du 



L . * Wiu-sn K . , L ' [^ 

/ — ( / , )du = — > / — = 

J** TTS n tl Jv /l - (sn/uR) 2 7rs„tl«^i ^1-sn/uR) 2 

= ^i:^(l + m + 2)-n 2 -^l(l + l)-n 2 ) 



n=—l 

First, in order to compare the cylindrical and spherical films, we take their areas to be the 



same; thus 2ttRL = 4ttR 2 or L = 2R. Then, if we define ^1(1 + 1) — n 2 as f n {l), a function 
of variable 1, the result of integration can be written as f n (l + l) — /„(/), which (for large 1) is 



the derivative f' n {l). In the limit of large 1, f n {l) ~ \Jl 2 — n 2 , whose derivative is l/y/l 2 — n 2 . 
Converting the sum to an integral the number of states becomes: 

2 1 1 
N - - V 

1 v states / j 



7T n= _ x V/ 2 - n 2 
dn 



71 J-l 



VI 2 - n 2 
= 21 

Therefore, at high frequency, the integrated number of states in the interval du is the 
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same in all three cases. 
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pictures 

1. The density of states of the cylindrical film (dashed curve) relative to that of planar 
film (dotted line) as a function of the reduced frequency ui R/s. The full vertical lines indicate 
the delta functions g sp h{^) for the spherical film. 

2. The heat capacity of a cylindrical (dashed) and spherical (full) film relative to that of 
a planar film, as a function of the reduced temperature defined in Eq. 13. 
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